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We elaborate on a powerful tidal interaction formalism where the multipole dynamics is kept 
generic and encoded in a linear response function. This response function is the gravitational coun- 
terpart of the atomic spectrum and can become of similar importance with the rise of gravitational 
wave astronomy. We find that the internal dynamics of nonrotating neutron stars admit a harmonic 
oscillator formulation yielding a simple interpretation of tides. A preliminary investigation of the 
black holes case is given. Our results fill the gap between Love numbers and dynamic tides. 
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I. INTRODUCTION 

Analytic models for gravitational interaction of com- 
pact objects in General Relativity (GR) are plagued by 
potentially very complicated internal dynamics. Recent 
progress on such tidal interactions is mostly focused on 
nondynamical models [1-3], which in particular can not 
account for oscillation modes. This situation was already 
criticized and improved in [4, 5]. However, it may be dif- 
ficult to extend this approach to objects other than neu- 
tron stars (NS) and the internal dynamics is developed 
around the Newtonian limit (but the adopted GR correc- 
tions seem to be sufficient for most applications). Here 
we devise a substantially more powerful tidal interaction 
formalism based on an effective field theory (EFT) ap- 
proach [6]. This approach was proposed in the context 
of black hole absorption [7] and consists in effectively 
replacing the extended object by a point particle com- 
prising dynamic covariant multipolar degrees of freedom 
(DOF). In this article the dynamics of the multipoles is 
kept generic and encoded by a linear response function 
to external tidal fields. 

Motion of extended bodies in General Relativity (GR) 
has been subject to question from the very beginning 
of the theory and gives the most important way to test 
gravity. Describing the dynamics of these objects is com- 
plicated and approximate methods have been developed, 
such as multipole expansion schemes along the lines of 
Mathisson, Papapctrou, and Dixon [8-11], between many 
others. However, the definition of covariant compact- 
source multipoles in GR according to Dixon is only useful 
for test bodies. The extension to self gravitating objects 
is not fully understood, though it is clear that Dixon's 
multipoles should be renormalizcd [12]. 

The adopted EFT approach implies a definition of co- 
variant source multipoles of self-gravitating objects in 
GR. This definition is implicit until an explicit matching 
of the point-particle description to the actual extended 
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object is worked out. This is the main purpose of the 
present work, but also highly nontrivial. (For instance, 
the determination of NS multipoles is already quite sub- 
tle for nonperturbed stationary spacetimes [13], see also 
[14].) This article substantially improves the situation 
for linear perturbations around a static (nonlinear) back- 
ground such as NS, Black Holes (BH), White Dwarfs and 
likely even Boson Stars. 

We illustrate our formalism with a simple Neutron Star 
(NS) model. The tidal constants (Love numbers and yet 
undetermined constants) are easily extracted from the 
response function. It turns out that, as long as linear 
perturbations are applicable, the internal dynamics of NS 
admits a formulation in terms of harmonic oscillator am- 
plitudes [15-17] similar to the Newtonian case, making 
the multipolar DOF composite. This astonishing result 
leads to simple and intuitive interpretations of tidal in- 
teraction in GR. A full analysis of the black hole case is 
still in progress. Its outcome is hard to foresee and thus 
for sure will bear surprises. 

Although numerical simulations capture the nonlin- 
ear aspects of tidal interactions, complementary analytic 
models stimulate invaluable (at least qualitative) inter- 
pretations of the physical processes at hand. An impor- 
tant aspect of our analytic dynamic tidal model is to 
naturally account for resonances between external tidal 
fields and oscillation modes of the NS in GR (see [16, 18- 
25]). Such resonances are of great importance. For in- 
stance, it was suggested recently that the oscillations ex- 
cited by these resonances can be strong enough to shatter 
the NS crust, thus producing a weak short Gamma Ray 
Burst (GRB) [26] (more precisely, a weak precursor to 
the main flare of the GRB produced by the merger of the 
binary). Besides such spectacular effects, resonances can 
of course leave more subtle, but invaluable, imprints on 
the internal structure in the Gravitational Wave (GW) 
signal. Numeric relativity simulations suggest that os- 
cillations excited by resonances can even be driven into 
the nonlinear regime and thus contribute signiflcantly to 
GW [27]. 

The next revolution in GR will certainly arise from GW 
observatories like Advanced LIGO and VIRGO. These 



detectors will begin its operation soon and likely detect 
GW from binary NS mergers on a regular basis [28] . Such 
GW signals encode a tremendous amount of information 
on the internal structure of NS. This expectation is sup- 
ported by recent numerical simulations, which reveal im- 
prints of the equations of state [29, 30] or the formation 
of a metastable hypermassive NS [31]. Simultaneous de- 
tection of GW and GRB can provide for the first time 
persuasive evidence for certain GRB scenarios [26, 32]. 

The present article is a continuation of our work in 
[17] (on the Newtonian case) and we adopt notations and 
conventions therein. 



II. EFFECTIVE ACTION 

Our approach to account for the innumerable internal 
DOF of compact objects follows along the lines of ther- 
modynamics. The obstacle is to identify state variables, 
which by definition describe the system on large scales 
(infrared, IR). In the case of gravitational interaction, 
these state variables reduce to the source multipole mo- 
ments: At the same time they encode the IR field and 
the motion [11] of the object. 

An effective point-particle action along the lines of [6, 
7, 33] is most natural to implement covariant multipoles 
as macroscopic variables. 
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where m is the mass of the NS and Eat is the electric part 
of the Weyl tensor. For simplicity, we only discuss the 
covariant electric type quadrupole Q"'' here, but inclu- 
sion of other multipoles is straightforward, see [33, Eq. 
(1)]. The indices a, b indicate the spatial components in a 
local Lorentz frame comoving with the NS. The worldline 
parameter r is the proper time here. 

We consider linear perturbations of compact objects, 
so we expect a linear response of the quadrupole to the 
(quadrupolar) tidal field Eab, 
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where the tilde denotes Fourier transformation from r to 
w, and F is the linear response function (or propagator). 
The main objective of the present article is to determine 
F from a matching procedure. As explained in [17], from 
a Taylor-expansion 
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the tidal constants /i2 and /ij emerge. The first pa- 
rameter /^2 is related to the dimensionless (relativis- 
tic, quadrupolar, 2nd-kind) tidal Love number fc2 = 
3GiJ.2/2R^, where R is the radius and G is the Newton 
constant, in agreement with definitions in [2]. Further- 
more, H2 parametrizes the tidal response beyond the adi- 
abatic case. Though it was formally introduced in [34], it 



was not determined numerically yet. It obviously comes 
out as a byproduct within our approach. The constant 
A is related to absorption [7], see also [35-37] and for 
a non-EFT treatment see, e.g., [38]. The time depen- 
dence of the mass parameter in the effective action due 
to absorption is discussed in [37]. The response func- 
tion is analogous to the refractive index in optics, where 
imaginary parts also encode absorption. This analogy 
enlightens the matching procedure. Indeed, the phase 
shift between ingoing and outgoing waves encodes the 
real part of the response (Love number/refractive index) 
while the change in amplitude is due to absorption. But 
the nonlinear nature of GR makes the interpretation of 
phase shifts more subtle. 

We should stress that besides encoding all quadrupolar 
tidal constants in a single function F, our approach can 
naturally accommodate the presence of oscillation modes 
that are obviously missed by a Taylor expansion (3). This 
possibility was not discussed in [7] , where the focus is on 
absorption. 

Generic extensions of the point-mass action were con- 
sidered in [39] and the resulting EOM were related to 
Dixon's results. This can readily be applied to (1). Ex- 
plicit expressions for the stress tensor in terms of Dirac 
delta distributions can be found in [40, 41]. However, it 
should be noted that the relation between Dixon's co- 
variant multipole moments and the covariant moments 
used in the action (1) is more of a formal nature when 
self-gravitating objects are considered. 

It is straightforward to derive the Newtonian interac- 
tion potential for binaries belonging to the effective ac- 
tion. Even the first post-Newtonian (FN) correction for a 
generic quadrupole was already worked out [42] (though 
not from the effective action; see also [43] for the impact 
on GW). However, the dynamics of the quadrupole was 
essentially left open and only made explicit for the adia- 
batic case. The present work fills this gap by providing a 
dynamical quadrupole model. It should be stressed that 
even if the effective action is applied to Newtonian or PN 
approximations, the response function F encodes strong 
field aspects of GR. This is the eminent advantage of the 
EFT approach. PN interaction potentials including tidal 
coefficients were derived in [34, 42, 44]. 



III. PERTURBED COMPACT OBJECTS 

Without going into detail, we just mention here that 
perturbations of static compact objects can be deter- 
mined from a system of coupled ordinary differential 
equations with the radial coordinate as the variable and 
the frequency a; entering as a parameter. (The specific 
case of nonrotating spherical symmetric NS perturbation 
goes back to [45, 46], for reviews see [47-49].) In the 
exterior, the perturbation equations are given by the fa- 
mous Zerihi [50] or Regge- Wheeler (RW) [51] equations 
for electric- or magnetic-parity type perturbations, re- 
spectively. These are the same equations that describe 



perturbations of Schwarzschild BH. The Zcrilli equation 
can be cast into the (simpler) RW form [52], such that the 
discussion can be restricted to the latter. Moreover, the 
RW equation possesses analytic series solutions [53], see 
also [54, 55], which are central for the present work. Our 
approach consists in solving the perturbation equations 
numerically in the interior and connecting to the analytic 
vacuum solutions by imposing appropriate boundary con- 
ditions at the surface. 

As the RW equation is a second order homogeneous 
differential equation, its generic solutions can be repre- 
sented by a linear combination of two independent so- 
lutions. In [53], the pair of analytic solutions Xq and 
Xq^~^ represented by series of Gauss Hypergeometric 
Functions converges at the BH horizon but not at spatial 
infinity, while for X'^ and X'^'^~'^ (series of Confluent Hy- 
pergeometric Functions) it is the other way around. Here 
v is the renormalized angular momentum, which is fixed 
by requiring convergence of the analytic series solutions 
[53] . We review the needed elements of [53] in Appendix 
A 2. The solutions can be matched as 
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where K^, is given by [53, Eq. (4.2)] or (A34). 

In the absence of dissipation, F should be real so it is 
natural to work with manifestly real quantities. Before 
proceeding, we thus introduce normalization constants 
TVjy such that the analytic solutions X'^ := N,^X^ are 
real. We also require that the asymptotic amplitude is 1, 
i.e., Ni, is uniquely defined by the requirement 
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where r^ = r + 2M\og{r/2M — 1) is the tortoise coor- 
dinate, r is the Schwarzschild radial coordinate, and we 
identify M = Gm. It is straightforward to work out an 
explicit series representation for N^, from the formulas 
provided in [53, 54]. The result is shown in Appendix 
A 2 c. Now the RW function X in the exterior can be 
decomposed as 



X = AiX!^ 
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The main numeric result needed for our investigation is 
encoded in the amplitudes Ai and A2. The introduction 
of e^ is suggested by an analysis of K^, and K-i^^i for 
small e. 

We checked our implementation of the analytic solu- 
tions [53] against a direct integration method [56]. 



IV. EFT CALCULATION 

The generic idea is to replace the compact object by 
an effective source encoded by the action (1), such that 
the RW function X at large distance (in the IR) is repro- 
duced. This singular source can be expressed in terms 



of Dirac delta distributions. In the effective theory we 
therefore need to solve an inhomogeneous RW equation 
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where I is the angular momentum quantum number. The 
source term S can be derived by projecting the stress 
tensor following from the action (1) onto tensor spherical 
harmonics, completely analogous to a point-mass source 
[57, 58]. Explicit expressions for the quadrupole case 
I = 2 are supplied in Appendix A3. The principle is the 
same for other values of I. 

As the distributional source S should mimic the com- 
pact object, it must be located at r = 0. However, due to 
the (regular) singularity at r = 2M, the inhomogeneous 
RW equation (7) then does not seem to make sense. This 
problem is resolved by understanding the solutions as 
expanded in the post-Minkowskian expansion parameter 
Af or e = 2AIuj. Expanding the solutions in e is subtle 
due to various poles arising from Gamma Functions. If 
one keeps I generic and performs the limit I — >■ 2 after 
the expansion, one ends up with a different set of solu- 
tions to the RW equation denoted Xg and X^''^^, where 
di = I — 2 represents the deviation from the quadrupole 
case 1 = 2. It holds 
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Notice that the 5,-poles in the coefficients are canceled by 
poles contained in the solutions X^ and X^'^^^ . Despite 

these complication we work with Xg and X^'^~^. The 
reason is that by keeping I generic one can easily identify 
the terms produced by the source S. 

Still the solutions may be singular at r = 0, so regular- 
ization techniques are needed to handle the delta distri- 
butions contained in the source S. We choose a smooth 
ultraviolet "cutoff' in the form of a Riesz-kernel repre- 
sentation, see, e.g., [59]. Furthermore, as the calculation 
of the source S for generic I is not an easy one, we use 
a more ad hoc approach. We take S for / = 2 only and 
multiply the Riesz kernel by r~ ' to augment it with a 
fractional multipole character. Finally we represent 5{r) 

by 
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where c? = 3 is the number of spatial dimensions. fiQ and 
c; are arbitrary parameters of unit inverse length for- 
mally introduced to make the expression dimensionally 



correct. The limit (5 ^- is understood to be taken in 
final expressions. 

We are going to construct the solution to the inhomo- 
geneous equation from the solution to the homogeneous 
one using the standard method of variation of parame- 
ters/constants. This allows us to reinterpret the analytic 
solutions X^ and X^'^~^, which converge in the IR, as 
the ones belonging to the EFT. So the ansatz reads 



X — CiXl + 6*2^ ■;, 



(11) 



where Ci and C2 are yet unknown functions of r. If Ci 
and C2 were constants, then this would just be a solution 
of the homogeneous equation. It is an elementary result 
that 
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where Hi and H2 are yet undetermined integration con- 
stants and W^, is the Wronskian w.r.t. the tortoise coor- 
dinate, 
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By virtue of Abel's identity, this Wronskian is actually 
a constant and can be evaluated at r = 00 by analyz- 
ing the asymptotic behavior of the analytic solutions, see 
(A53). This procedure to solve the RW equation with 
a delta source has some similarities to the construction 
of the gravitational Green function, which was also ob- 
tained from analytic solutions very recently [60]. It can 
be interesting to further study this connection in the fu- 
ture. 

Next one must constrain the integration constants Hi 
and H2. The solution they represent is not allowed to 
correspond to further delta-type sources at r = 0. This 
argument leads to H2 = 0, as the sourced terms are pro- 
portional to r~'~^ ~ Xg'^~^, while Xg ^ rK If we would 
set I = 2 in the very beginning, then one must identify 
the constraint on Hi and H2 by iteratively solving the 
RW equation in e. However, the coefficients in the nu- 
meric solution (6) scale differently by 4 orders in e, so one 
must iterate to the same order in e just for the leading 
order result. This is the reason for keeping I generic here, 
at the cost of introducing poles in 5i into the calculation. 
Dimensional regularization would solve this problem in a 
similar manner, see [61]. Interestingly in the static limit 
the dependence on the dimension can be absorbed into I 

[61]. 

It should be emphasized that the lower integration 
bound in (12) is r = 0, which is only possible if the inte- 
grand is understood as expanded in e. Then the singular- 
ity of dr^/dr and the oscillatory behavior of the analytic 
solutions at the horizon are removed. Here nonconver- 
gence of Xg and X^''~^ at the horizon is actually not a 



flaw, but a feature. Furthermore, Eqs. (12) only start to 
depend on the upper integration bound r' at linear order 
in S, so Ci and C2 are actually constant for a; ^- 0. 

Finally, the quadrupole components entering the 
source 5" must be determined according to (2). This re- 
quires to evaluate the tidal field E°-^ at r = 0, or 
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where again the Riesz kernel provides the necessary regu- 
larization and the integrand is expanded in e before eval- 
uation. Note that E°-''{r,Lu) follows from the RW master 
function (11) and the background. See Appendix A3d 
for its representation in RW gauge. 



V. RESULTS 

Comparing the numerically obtained solution (6) with 
the effective one (11) leads to two conditions. These are 
solved for the yet undetermined quantities Hi and F in 
terms of Ai and A2, completing the computation. At the 
end of the day, one arrives at 
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Let us insist again that though this seems to be an expan- 
sion in e, the numerical quantity A2/A1 can still capture 
strong field effects from the interior. We define a renor- 
malized Fms by dropping the poles in 6 and Si , analogous 
to minimal subtraction (MS) in dimensional regulariza- 
tion. 

We apply our method to the astrophysically most rel- 
evant case of NS. The used system of differential equa- 
tions is derived in Appendix A 1. We use a simple poly- 
tropic equation of state (EOS) with index 1 for the nu- 
clear matter. The results presented below are for a NS 
with m = 1.277710 and R ~ 8.85 km. The complex quasi- 
normal mode frequencies for this specific NS model were 
reported in [47] and our numeric implementation repro- 
duces them very well (except for the damping of the cur- 
vature modes). 

An excellent fit for Fms for this NS model turns out 
to be 
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provided we also fit the renormalization scale ^o- We nu- 
merically generated a set of 350 data points (with higher 
density near the poles) and the fit deviates from all of 
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FIG. 1. Propagator of the quadrupole in minimal subtraction 
(MS). The dots are some selected numeric data points, the 
solid line is the fit (16) and (17). For the circles the loga- 
rithmic scale dependence was ignored (/io ~ oj) and for the 
crosses only the leading order in (15) was taken into account. 
This shows that the e^ corrections in (15) are essential for 
a good fit, while the corrections from /io start to contribute 
only beyond the f-mode. 



them by at most 2% (sec Fig. 1). The optimal fit param- 
eters are given by 

w/i? = 27r 0.0851, (7/ = 1.98 X 10"^ (17) 

Wpi? = 27r 0.194, Qp = 9.1 X 10"'*, ^lQR^Q.&. 

It is straightforward to infer the tidal constants defined 
by (3). It should be noted that all fit parameters are 
essentially independent of c; , which we varied from uj to 

It is remarkable that the quadrupole propagator can 
be approximated by a sum of response functions of har- 
monic oscillators, just like in the Newtonian case [17]. 
The relativistic case thus seems to admit an amplitude 
formulation [15, 16] analogous to the Newtonian case. 
The a priori very complex internal dynamics of the NS 
is then approximated just by a set of harmonic oscil- 
lators, which are the more fundamental effective DOF 
composing the dynamical quadrupole. The constants 
qfB? /G and QpR^/G can be understood as GR versions 
of the overlap integrals. Resonances are quantitatively 
described by forced harmonic oscillators. An extension 
of this mechanical picture to nonlinear oscillators and/or 
mode coupling can offer a demonstrative phenomenolog- 
ical way to model even more realistic situations. 

The frequency dependent Love number introduced re- 
cently in [62], where the formalism in [4, 5] was used, 
should be related to our response function in the low fre- 
quency regime and a comparison is most interesting due 
to the very different setup (single object in perturbation 
theory vs. complete binary). It is further an interesting 
question whether the expansion (3) including the tidal 
coefficient ^'2 introduced in [34] is enough to find agree- 
ment within the regime where a comparison is possible. 

It is well known that the f-mode properties (frequency) 



basically just depend on the mean density of the NS [47]. 
To a good approximation this should also be true for the 
tidal response of the quadrupole, as we find here that 
it is largely dominated by the f-mode (at least for the 
adopted NS model). This aspect may be related to the 
universal relations discussed in [62-64]. 

Also the BH case can be readily investigated using (15). 
One can even work out analytic formulas for Ai and A2 
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[53], reading 
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The nontrivial character of this analytic result becomes 
apparent once (A34), (A48), (A46), and (A47) are in- 
serted. Because of the absorption due to the horizon 
one can not expect poles for cj e M like in the NS case. 
Tliis makes the analysis more complicated, as the whole 
complex plane must be considered. However, one can 
immediately obtain an expansion of _F in e, 
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This implies that the BH Love number /X2 vanishes, in 
agreement with the findings in [2, 3, 61, 65]. Unfortu- 
nately the unspecified parameter ci can substantially in- 
fluence the next order tidal coefficient /ij. This makes 
clear that for a rigorous investigation one should first 
redo the EFT calculation within a better regularization 
method like dimensional regularization. 

But it should be stressed that the Love number is inde- 
pendent of ci. In fact, it is possible to obtain the leading 
order of (15) by setting Si = throughout the compu- 
tation using a shortcut. Though I is not available to 
identify the constraint on Hi and H2 , a simple argument 
can be given at leading order. Hi and H2 must corre- 
spond to the linear combination of X^ and X^'^~^ that 
cancels the z~^-term in the solution at orders e^ and e^. 
Unfortunately at order e^ [corresponding to e^ in (15)] 
this approach breaks down due to z~'^ log z contributions 
to the solution, which can only be interpreted by an iter- 
ation of the field equations. Still it is highly desirable to 
reach the next to leading order, as illustrated by Fig. 1. 



VI. CONCLUSIONS 

The validity of our results is supported in many ways. 
First of all, by comparing the adiabatic limit of the 
analytic solutions used here against the analytic zero- 
frequency solutions used in [2] one can show that the 
definition of //2 in [2] agrees with the definition through 
(3) and (15). Second, the term linear in e in (19) was al- 
ready derived in [7] and agrees with our findings. Finally, 
an intermediate result is the dependence of quadrupole 
components Q"-^ from (14) on the scale fiQ, which in fact 



agrees with the beta function found in [33] using dimen- 
sional regularization (if it is assumed that c/ is indepen- 
dent of /Xo ) . 

But unfortunately our current results depend on an- 
other parameter c; with unclear interpretation. This sit- 
uation can be improved by applying dimensional regular- 
ization to the EFT calculation, which is also most use- 
ful for applications to post-Newtonian theory. But we 
expect our results to be good approximations for the di- 
mensional regularized ones (the matching scale /ip will 
be slightly different). Still this calls for a clearer connec- 
tion between our formalism and Ref. [33] (where also the 
background Schwarzschild metric is generated within the 
EFT). At the same time, higher multipoles (including 
magnetic/axial) should be treated as well. 

However, these current problems with the regulariza- 
tion method play no role for the static limit. Thus no am- 
biguities for the definition of the Love numbers emerge. 
The predictions for the RW function from the EFT are 
"simply" matched to the numeric results, as in [61] for 
the black hole case in dimensional regularization. Conse- 
quently there is no need to interpret the definition of fc2 as 
relative to the BH case in the current approach. Further- 
more, our computation is based on the exterior solution, 
so it is applicable to arbitrary (nonrotating) compact ob- 
jects. 

Another obvious next step is an application of our 
method to more realistic NS models. Besides realistic 
EOS, an investigation of the NS crust is most promising 
due to a possible connection to precursors flares in short 
GRB [26] by a shattering of the crust. 

Further, realistic NS are rotating. Neutron star modes 
can become unstable in the rotating case (including the 
f-mode [66]), which can also give rise to violent astrophys- 
ical processes. An extension of the present method can 
be tried within a slow rotation approximation, see, e.g., 
[67, 68] and references therein. The r-modes of rotating 
NS are of particular interest for resonances [16]. 

Analytic predictions for GW including tidal effects 
from an Effective One-Body (EOB) approach agree even 
quantitatively with numeric simulations [44, 69] . Yet the 
tides are modeled by Love numbers and absorption [70] 
only. Inclusion of the dynamical multipole response func- 
tion into the EOB formalism is expected to establish the 
impact of resonances between NS modes and orbital mo- 
tion on GW in a reliable manner. 

Our approach can probably also be evolved into a 
method for finding oscillation modes. In the conserva- 
tive case, modes can be found by "just" integrating the 



perturbation equations for real frequencies. Estimates for 
the damping times can be obtained using the quadrupole 
formula (and basically correspond to the overlap inte- 
grals). However, if the modes are damped by, e.g., dis- 
sipative effects in the nuclear matter or mode coupling, 
then the poles of the response function should have a 
nonvanishing imaginary part. This also illustrates that 
our approach separates properties of the star from effects 
due to the surrounding spacetime, which is nontrivial as 
the background is nonlinear and does not admit super- 
position arguments. 

The BH case is largely left unexplored for now. At 
the same time, the prospects are fascinating. If a me- 
chanical oscillator model for BH can be formulated, then 
one can further elaborate on the thermodynamic analogy 
of these macroscopic DOE. This can lead to insights on 
macroscopic concepts like BH entropy and temperature 
from an EFT point of view. Analogies to the AdS/CFT 
correspondence discussed in [7, 71] can probably be made 
more explicit, too. The highly damped modes and branch 
cuts (eventually introduced by nonanalytic terms like 
loge) can be difficult to handle. Extension to the case 
of rotating BH should be almost straightforward, as the 
perturbation master equations [72] , the appropriate ana- 
lytic solution [73] , and the effective action [35] are readily 
available. The Love number of rotating BH should even 
come out unambiguously if the current (improvable) reg- 
ularization method is applied. 

For an application to BH scattering a simplistic fit 
of the response function can be accurate enough. This 
can lead to interesting connections to, e.g., the scat- 
tering thresholds discussed in [74, 75]. Finally, if the 
present method is applied to perturbations of massive 
scalar fields around rotating black holes, the fioating or- 
bits existing for extreme mass ratios [76] can possibly 
be constructed for comparable mass binaries using PN 
methods. 
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Appendix A: Formulas and Implementations 
1. NS Perturbation Equations 

a. Preliminaries 

A static spherically symmetric star configuration is described by a metric 

dsl = -f{r)dt^ + b{r)dr^ + r^d^^, (Al) 

where dVt is the line element on the unit sphere, and a perfect fluid with stress tensor is given by 

T^, = {p + P)u^u, + Pg^,, (A2) 

where p is the density, P is the pressure, g^^i, is the metric components, u^ is the four- velocity such that u^w^ = — 1, 
and Greek indices belong to the spacetime coordinate basis. In the NS case, the description is usually completed by 
a barotropic EOS relating P and p (neglecting temperature), according to 

p(r)=p(P(r)), (A3) 

for a given function p. In this article we considered a polytropic EOS defined by 

P^Kp"^. (A4) 

The dynamics of the compact object is then given by the Einstein equations and the conservation equation of the 
stress tensor 

Rf.. - \R9f.u = IGttGT^,, V^T^" = 0. (A5) 

The even parity metric perturbations around a spherically symmetric background are given in the Regge- Wheeler 
gauge by 

dsl ^ ~f{r)ho{x^')dt^ + 2iLuhl{x^')drdt + 6(r)/i2(.T^)dr2 + r'^k{x'')dn^, (A6) 

where the functions h e {/ig, /ii, /12, k} depend on the coordinates x^ according to h = h{r) cxp{iujt)Yimi^) , Yim being 
the scalar spherical harmonics. The total metric then reads ds^ = ds^ + dsf, where the subscripts and 1 denote the 
background and perturbation, respectively. The perturbations to the matter fields are given by 

P = P„{r) + P,e-''^*Yi„,{n), p = p„{r)+p,{r)e-''^'Yi^{n), u^ = <(r) + <(r, 0, V5)e~'"*. (A7) 

Given an equation of state, the perturbation to the density is given in terms of the perturbation to the pressure: 



dPo ' c2 



PI = 377^1 = ^. (A8) 



which defines the speed of sound Cs ■ 

In the following, it will be useful to introduce the function U such that 

Pi=(u-^yp, + Po). (A9) 

The solution for u'^ is given in terms of the function U and of the metric perturbation 

- Vmu^^dx, = (^dt + ^d^in^immdr) v^n) + W^v^Y.^dn., (aio 

\ 2 ujb[r) J r-^oj 

where dVti — [dO, dip). 
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b. Master Equations 

The perturbed equations can be solved for Hq, /i2 and their derivatives, leaving three ordinary differential equations 
in the interior of the compact object 

U" = -{b\fk{r + 8nr^ ^)'(-l + p) + Qrf{f U'{1 + p') + hi{b + p')) - b {f{2f{W{l{l + 
-167rr2(P + p)) + r(-8 + I + l^ + 8Trr'^{-7P + p))U') - ir^k{^l + p) + 2r(-6r uj'^U + f{A + l + f 
+32TTr^P)U') p') + rhi{-Ar^uj'^{-l + p') + /(16 + 11/(1 + l) + {8 + l + f)p' + 327rr2 (P - 3p + 2Pp')))) 

+ (6 f{Ar'^Lo^k{-l + p') + 2 (/)^(2C/(1 + l + f- %TTr^p){l{\ + l)~ IQ-Kr^iP + p)) + rt/'(-3(l + Z + ^2) 

+8 7rr2(-P(3 + 2/(l + /) + 87rr2p) + (3 + / + ?2+ ^-kt^ P)p - ^-Kr"^ {pf) + p 

+ {l + l^ + 87rr2p(2 + / + /^ + 87rr2 P))p')) + r/(/ii(l + 87rr2p)(2 + 31 {1 + I) + 167Tr^{~2p + P(-l + p')) 

+ (2 + 1 + P)p') + 2r{-k{l + l' + Wirr^ P)(-l + p'} + 2oj^U{8TTr^{P + p) ~ {1 + I + P + Sirr^P) p'))))} 

/{4rV'(3 -b{l + l + f- Sirr^p))}, (All) 

k' = -{327rrw26[/(p + p)-\k{-6f + {bff{l + 87rr2p)(-/(l + l) + Snr^P + Snr^p) 

+b{Ar'^uj'^ + /(6 + / + ^2 - 87rr2p - 2'iTTr'^ p)))] / [r f] - [hl{Ar'^uj'^ + /(/ + f - IGnr^P - 16T:r'^p){~2 

+b{l + P - Snr^P - Sirr'^p)))] /[r^f] - 167r/(P + p) (2 + 6(-Z(l + /) + Sirr^P 

+87Tr^p))U'}/{2{3 -b{l + l + l^~ Sirr^p)) } (A12) 

h[ = {6//1I - r{bffk{l + 8nr^pf - (5)^(4r^CJ^fc - /(2r(167rr VC/(P + p) + fc(l + 21 + 21^ + IQiir^P - Sirr^p)) 

+hl{2 + 31 + 3f - 128n^r^{Pf - 8{5 + I + f)TTr^p + QAir^r^ipf - 87rr2p(l -21- 2f + 24:nr^p))) 

+167rr(/)2(P + p){2U{l + 1 + P - Sirr^p) + r(l + 87rr^P)U')) + b{-hl{Ar^u^ + /(8 + 3/ + 3f 

+87rr2p - 56TTr^p)) + rf{-9k + 167r/(P + p){6U + rU')))}/{2rf{3 - b{l + I + P - STrr^p))}, (A13) 

where we omitted the radial dependence of the functions. 

These equations can be further written in terms of the Rcgge- Wheeler Master function X by using the following 
change of function 

k = {{l{l + 1){-24:M^ + l2Mr + {-I + 1)1(1 + l){2 + I) r^) -24:MrW)X 

+2l{l + l)r{-2M +r){6M + {-2 + l + P)r)X'}/{2l{l+ l){-2 + I + P)r^}, (A14) 

hi = {(-72A/3 + 12(3 + I + f)M2r - (-1 + 0^(1 + 0(2 + 0^^ + 3Mr^{l{l + l){-4 + I + f) + 4r^uj^))X 

-{2M ~ r)r{6M - l{l + l)r) (6M + (-2 + / + l^)r)X'}/ {l{l + l){-2 + l + f)r^{-2M + r)}. (A15) 

This transformation has the advantage to turn the vacuum equation {U — 0) to the simple Rcgge- Wheeler equation 
(in tortoise coordinates) 



where r, is the tortoise coordinate defined by 



^ + ''^^°S(2^-^)- ^^''^ 



Note that the Reggc- Wheeler equation originally describes the odd sector of the metric perturbations. The system 
of equations describing the even sector is called the Zerilli equation. However, these sectors are isospectral and are 
actually equivalent, by virtue of the transformation (A14). Here we use the Regge- Wheeler equation since the form 
of it is of Heun's equation and is better suited for the construction of the analytic series solutions of [53] . 

c. Boundary Conditions and Series 

Within the NS interior it is difficult to extend the numeric integration up to the boundary points r = and 
r — R. This problem is solved by terminating the numeric integration very close to these points and use analytic 
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series solutions to extend the numeric solutions to the boundaries. Furthermore, certain boundary conditions must 
be fulfilled, which are directly implemented into the series solutions here. 

At the origin r — 0, the relevant boundary condition is simply the regularity of the perturbation master functions. 
This imposes two independent conditions, so we need only two (of four) integration constants to parametrize the 
solution around r = 0. The leading order series solutions read 

U{r)^Uor'[l + 0{r)], (A18) 

k{r) = kor'[l + 0{r)], (A19) 

hi{r) = -^(fco - 87rG;7o/(0)[P(0) + p(P(0))] + 0(r)), (A20) 

where Uq and fco are the integration constants. 

The boundary condition at the surface is given by the requirement that the Lagrangian (comoving) perturbation 
of the pressure or density vanishes for r = i?, or explicitly 

U'\n = -fk{3f + b'fP' + b{4RW - 2/((/ + 2)(/ - 1) + 2P))) + j^0L 

+ ^ (-6/ - b^fP{2p + l{l + l))+ b{-AR^oj^ + f{8P + /(I + 0)) ) 1^ , 

D = 2f{bP-l){biP + l{l + l))-3) , P=l + 8nGR'^P. (A21) 

However, the behavior of the perturbations near the surface r = R crucially depends on the EOS. In the following 
analysis, we restrict to the case that the EOS near the surface is a polytrope with index 1 < n < oo. Then the 
boundary condition is actually equivalent to regularity of the perturbation master functions. The boundary condition 
allows us to eliminate one of the four integration constants, so we are left with U{R), k{R), and hi{R). At the surface 
the functions k and hi must be continuous, which provides two further boundary conditions. We are therefore able 
to express k{R) and hi{R) in terms of the RW function X{R) and its derivative X'{R) (which describe the exterior 
perturbation). The boundary series finally reads 

U(r) = U„ 

+ (-1 + /);(! + 0(2 + l)R'^{l + 1^ - 2R^iu^) + 2AM^R^{1 + f + R^lo^) 
- 2MR^{{-1 + l)l{l + 0(2 + /)(1 + l + r^) + 12i?2a;2)) 

+ MR^{-{-l + l)l{l + 0(2 + /) - 12i?V)) I 

2 f (-1(1 + l)M{2M ~R) + {M- 2R)R^uj^)Ur 



+ {R-r) 



MR{-2M + Rf 
Xr 



(5761'/^ - 96(3 + I + f)M^R 



2l{l + /)(-2 + I + P)M{2M - RfR^ 
- 3(-l + 0^(1 + 0(2 + /)(2 + / + f)MR^ + 2(-12 + (-1 + l)l{l + /)(2 + l))MR^oj^- 
+ (-1 + 0^(1 + 0(2 + 0^'^(^ + '^ - 2i?,2cj2) - 12M^R^{1{1 + 0(-10 + 3^(1 + 0) + 14i?V) 
+ 2Af2i?3((_i + ^)/(i + /)(2 + 0(13 + ? + f) + 48i?V)) 

VI 

^ :(144M^-48M3i?- (-1 + 0^(1 + 0(2 + 0^" 



l{l + l){-2 + l + P)MR{-2M + Ry 
+ Af2i?2(_7(_i + /);(! + 0(2 + - 12i?^a;2) + AMR^{{-1 + l)l{l + l){2 + I) + SR'^lo^)) | 
0[(i?-r)3], (A22) 



k{r) = 
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(/(I + 0(-24A/2 + 12MR + (-1 + 0^(1 + 0(2 + l)R^) - 2AMR^uj^)Xii (2M - R){QM +{-2 + 1 + P)R)X'ji 
2l{l + l){-2 + l + P)R^ {-2 + l + P)R^ 

+ ^^ "^'[K 2i?2 + (-2 + / + /2)(_2M + i?)J^^+V 2R +/(l + /)(-2 + Z + ;2)j^«j 

+ O [(i? - r)2] , (A23) 

'^^^^ - lil + l)i-2 + l + P)RH-2M + R) ^-™' + 12(^ + ^ + ^'^'''^ 

- (-1 + /)/(! + 0(2 + 0^''^ + 31/i?2(/(l + /)(-4 + Z + /2) + 4i?2cj2)] 
(6A/ - l{l + l)R) (6il/ + (-2 + ; + P)R) X'j^ 

^ l{l + l){-2 + l + P)R 

+ (^ - < iil + l)i-2 + l + P)R^i-2M + Rr ^-'''''' + ''^' + ^^(^ + ^»^^^ 

+ (-1 + 0^(1 + 0(2 + l)R'^{l + 1^ - R^^^) + 12M^R^{1{1 + 0(-5 + l + f) + TR^uj"^) 

- MR^{1{1 + /)(-16 + /(I + 0(1 + 2^(1 + 0)) + 24i?2a;2)) 

X' 

^ ;{-72M^ + 12(2 + I + P)M'^R 



l{l + l){-2 + I + P){2M - R)R'^ 
- 2(-l + 0^(1 + 0(2 + 0^^ + MR^(1{1 + l){-20 + 7/(1 + 0) + 12i?2cj2))j 
+ 0[{R- rf] , (A24) 

where Uii,Xji,X'j^ are the functions U,X,X' evahiated at the radius R. 

If the RW equation is solved numerically using a direct integration method [56] , then one must derive series solutions 
for the RW function X at r = (X) in a similar manner (and for BH also at the horizon, where the physical boundary 
condition only permits an ingoing flux). 

2. Analytic Solutions from [53] 

a. Solutions to the RW Equation 

The pair of independent (UV or "near-zone") solutions {Xq,Xq'^^^} and the pair of independent (IR or "far-zone") 
solutions {Xp,Xj''~^} are given by [53, Eqs. (2.16) and (3.6)], see also [54, 55], 

Ao-e (-X) (1-x) 2.(1 X) a„ r(,i - ze + i. + 3) 

n— — oo 

X 2Fi{-n -ie-v-2,-n-ie~v + 2; -2n - 2v; 1/(1 - x)) (A25) 

X'^ = (l- IV" V 2-^+";"^-" -■'^+"+ia'^ T{n~i€ + v- l)T{n -ie + u + l) 
^ \ z) ^^^^ '^ T{2{n + v) + 2)T{n + ie + v + 'i) 

X iFi{n + ie + V + l]2{n + v) + 2;2iz) (A26) 



where 



e = 2Muj, (A27) 

z = wr, (A28) 

x^l^-^-1-f, (A29) 

.F,(a,6;c;z) = f i%i^^, (A30) 



n=0 



(c)„ n\ ' 



oo / ^ ^j 

iFi(a;c;z) = ^ffii^, (A31) 

n— 
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(a). 






(A32) 



Here r is the Gamma Function, 2F1 is the Gauss Hypcrgcomctric Function, iFi is the Confluent Hypergeometric 
Function, and (a)„ is the Pochhammer function (rising factorial). The coefficients a{^ entering the series and the 
renormalized angular momentum ly are explained in the next section. 
The relation between the solutions is given by [53, Eq. (4.1)] 



Xq = K^Xc 



Cj 



(A33) 



where [53, Eq. (4.2)] 



K,. 



- csc{tt{i' + ie)) 



E' 



r(r + ie + 1/ - l)r(r + ie + i^ + l)T{r + ie + ly + 3) 

r(n -ie + iy - l)T{n - ie + u + 1) 



, T{n + r + 2iy+ l)T{n + ie + u - I) 
'' {n-r)\T{n-ie + v + i) 

-1 



E <^ 



n— — 00 



(A34) 



' (r - n)\T{n + r + 2u + 2)T{n + ie + i/ + l)T{n + ie + i^ + 3) 
The value of r G Z is in principle arbitrary, which can also be checked numerically. For definiteness, we chose r = 0. 



b. Recurrence Relation for a'^ 

The three-term recurrence relation for the a^-^ reads [53, Eq. (2.5)] 

ana„+i + /^„a„ + 7„a„_i = 0, 
where the coefficients are given by [53, Eq. (2.6), (2.7), and (2.8)] 

^ ie(y + n — ie — l)(i^ + n — ie + l)(i^ + n + ie — 1) 



1: 



{v + n+l){2{iy + n) + 3) 

(e^+4), 



(v + n){iy + n + 1) 
ie{i' + n — ie + 2){i' + n + ie){v + n + ie + 2) 



2e' 



(l^ + 7l)(2(z/ + 7l)-l) 

We proceed along the lines of [53, 55] by defining continued fractions Rn(y) and Ln(y) [53, Eq. (2.9) and (2.10)] 

^n In 



(A35) 

(A36) 
(A37) 
(A38) 



R^{u) 



Ln{v) 






/3- + <fi«+i(i^)' 



From these expressions it is straightforward to infer that 



lim nRn{v) = -77^ 

n— >oo Z 



lim nLn{v) 



2' 



(A39) 
(A40) 

(A41) 
(A42) 



provided that the continued fractions converge in the specified limit. The corresponding solution to the three-term 
recurrence relation is called the minimal solution in the specific limit and is guaranteed to exist. But the minimal 
solutions for n — ?> 00 and n — > — cjo are not necessarily the same, e.g., in general one can fulfill either (A41) or (A42), 
but not both at the same time. However, requiring both (A41) and (A42) fixes the renormalized angular momentum 
V. This is dictated by the convergence of the analytic solutions to the RW equation. 

In practice, one uses the limit (A41) as a starting value for Rn{v) at some large but finite n > 0. From the continued 
fraction (A39) one can then easily determine Rn{v) for any smaller n. An analogous process can be applied to Ln{v), 
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this time starting from a large but finite negative n < 0. Finally, one imposes the consistency condition [53, Eq. 
(2.11)] 

i?„ML„_i(zy) = l, (A43) 

at some value for n. For definiteness, we chose n = 1. We solve this condition for i/ using standard numerical 
root-finding procedures starting from the initial value [53, Eq. (6.3)] 



p = l+ - 



(/-2)2(/ + 2)2 4 (?-l)2(/ + 3)2 \ £2 



2Z(2;- i)(2/ + i) /(; + i) (2Z + i)(2; + 2)(2; + 3) 



-^ 27TT + ^^^')- ^^^^) 



Finally, one can determine aj^, which is fixed up to an overall factor. As in [53], we set Oq = 1. Besides numeric 
approaches, it is of course possible to work out analytic series expansions in e for i> and aj'j, see [53, Sec. 6] for further 
discussions. 

c. Normalization of Analytic Solutions and Wronskian 

It is straightforward to determine the asymptotic behavior of X'^ as 

X'^ '•*^°° AJ^ine-*'^'^* + ^&oute'"''% (A45) 

with the complex amplitudes 

Aj;,„ = iz-+-i E -z"(2.)-e-(-+t)a- [:i"'" + "~^i , (A46) 

n— — oo ^ ^ 

^cont-^j Z. ^2e) a„^ ^^^^^^^ ^^^^^^ (A47) 

n— — OO ^ ^ ^ ^ 

Then wc obtain for the norniahzation 

N, = \{A'^c,^A^aout)-'^- (A48) 

Notice that N^ is multivalued (bi- valued). Eventually the second root must be used in some frequency regimes in 
order to make results continuous. In the present investigation this is necessary at about LoR/iir > 0.23. It follows 
that 

X^^ := N,X^, (A49) 



2 



^Cout \ -iujr, , I ^C out \ ibjr. 



^-i(u;r.+Q„) , i{u>r,+au) 



witha. :=:^log^^, (A51) 

^* ^C in 



1 

2 

cos(a;r* + a^), (A52) 



as envisaged. Because the RW equation has real coefficients, it is guaranteed that a real solution for a^ exists. 
It is also straightforward to obtain the Wronskian (13) from this analysis as 

W. = 2zuN,N^,^M'binAc7nl - ^&out^ch7'): (A53) 

= uishi(ay — a-y-i) (A54) 

The Wronskian (13) based on the solutions X^ and X'^^~'^ follows by expanding (A53) in e for generic I. 

3. Effective Source 

a. Inhomogeneous RW Equation 

The homogeneous Zerilli and Regge- Wheeler equations describe vacuum perturbations. Considering additional 
matter fields will source the vacuum perturbation equations. If the right hand side of the perturbed Einstein equation 
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is ST^i^, the combination of the metric perturbations leading to the source to Zerilli equation and further transformed 
to Regge- Wheeler equation is 

o ill 

^_ 2£2r2(-2M + r)'(6(-2 + £2)M+(4-2£2 + (-2 + 4)5£|((-2 + £2K2)2)r) 

— '-' — 9 Q Too 

-864M4 - 144(-5 + ^2)^^ + 36 (-4 + ef)A'Pr^ + 12 (-2 + £2) i2Mr^ + (-2 + £2) i2r^ 

1 1 1 o 

2V2£2(-4 + 2^2 - (-2 + £2) 2^1 ((-2 + ^2)^2)2 )r2(-2M + r)' 

/oi 



Toe 



(6A/ + (-2 + ^2)r) (-144A/3 + 72APr + 6(-2 + 4)^2Mr2 + (-2 + £2)2^27-3)0; 

4^(-2 + £2)^ ((-2 + £2)^2) ^r(-2A/ + r)'(-12il/2 + 12Mr + (-2 + £2)^) 

(6A/ + (-2 + e2)r) (-144Af3 + 72APr + 6(-2 + 4)£2Ai"r2 + (-2 + £2)2^^3)^ 

1 1 1 
2£2r^(6(-2 + £2)M + (4 - 2^2 + (-2 + £2) 2£| ((-2 + £2)^) 2 )r) 

-864A/4 _ 144(-5 + £2)]VPr + 36(-4 + £2)m2^2 ^ i2(-2 + £2)^^2Mr3 + (-2 + €2)^4^4 

2V2(-2 + £2) 2 ((-2 + £2)4) ^ (2Af - r)r4 ^ 

-144Af3 + 72Af 2^ _^ 5^_2 + £^)£2MP + (-2 + £2)^^2r3 

2V2(-2 + £2)M(-2 + ^2)^2)^(2A//~r)'r2 
(l44Af 3 _ 72il/2r - 6(-2 + £2)^2Af r2 - (-2 + 4)24r3)w °'^ 

2V2((-2 + £2)^2) ^ (2M - r)r4(6Af + (-2 + £2)r) ^ 

2 — r'*^' (Aooj 



-144Af3 + 72Af2r + 6(-2 + £2)£2MP + (-2 + £2)^2r3 

where we introduced the notation £2 ^l{l + l) and where Tzt Z & {00, 01, 11, Oe, Oo, le, lo, t, e, 0} are the ZeriUi tensor 
spherical harmonic (TSH) components of T^'' defined by 

Tz=Nz I T^'^Yl^^dn. (A56) 

Here Yz^fiu are the Zerilli TSH [57] and Nz their normalizations given by A^ = {1, —1, 1, —1, —1, 1, 1, 1, 1, 1}. 
Finally, the sourced Zerilli equation converted to Regge- Wheeler form is given by 



b. Stress Tensor 
The stress tensor up to the quadrupole approximation reads [40, 41] 

V^T^" = Jdr u^^p''^5(i) + ii?„0^(^ J'^'^^"(5(4) - V„(5"(^7/'')%)) - ^V0V„(J^(^'')"<5(4)) , (A58) 

where 

pM ^ ,^yM _ y^ + -U6i?ede'^ J'^^''^ (A59) 

as 6 

and (5(4) = 5{x^^ ~ z^). Here the 4-quadrupole J"''^'^ jj^s the same symmetries as the Riemann tensor 

jaPfif _ j[al3][fj,i^] _ jfifafi fA60) 

jMm]^ = ^ .r'P"" + J'^^"""" + J"""^" = 0. (A61) 

It results directly from the effective Lagrangian as [39] 

r^^- = -6^^, (A62) 



which is defined by (1), 



'off 



dr L\, 



L;, 



u /~\ah 



-TO- -E^^Ca^eti^Q 
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(A63) 



Here Ca'' is the tetrad defining the local frame. We formally extend the local spatial indices a, 6 by a time component 
here. This is fine if we also set all time components of quantities defined in the local frame to zero, e.g., Q"-^'^'^ = 0. 
For the sake of the variation, we can then consider Ca'^ as unconstrained. (We implement the constraint e'-^^'' = u^ at 
the level of the equations of motion.) Notice that the result (A62) from [39] is valid in the presence of a tetrad Ca^ . 
Using E^i, = R^avpu^u'^ (in vacuum) we get 

ja/3A.i.^_3y[ag/3][^^H^ (A64) 

which is what we anticipated. The spin vanishes here, S^^'^ = 0. We further disregard the mass term to, as we are 
only interested in the contributions from the quadrupole here. 



c. Quadrupole Source for RW Equation 



For technical reasons we are not working with a local Cartesian basis, but in one that is adapted to TSH, i.e., 

VTsn = .9p.e"''e''^ (A65) 

r/^^SH = diag(-l, 1,1,1/ sin^^), (A66) 

77jSH = diag(-l,l,l,sin2 0). (A67) 

We can then transform components in this local basis to TSH components in the usual way. Our choice for the frame 
field reads 

/ hajr) 



I — 7^ \ 



(e'^") 







/h(r) 











+ e-*"*y,„(r!) 



V —)-^ j 



\ 



^h('-) /i2(r) Q g 





V 



^ 



y^^/ 



Obviously it fulfills e(°)^ = u'". 

The required components of the stress tensor in TSH basis 7x are given by 

\/3Qic(r(r-2A/),5;(r) + (2r-3A/)<5,(r)) , ^nQMr) 



Tor 



r2(2M-r) ^r'^{r-2M) 

St (r ((5r - 7A/)5;(r) + r{r - 2M)5';{r)) + i{M + 2r)5r{r)) 






_ Qtio^/r - 2M{r5l{r) + 3(5,. (r)) V^Sio^'^'-lO 

Toi — T7T\ 1- 



r3/2 



rii = 



Toe 



QtVr - 2M {5r{r) {-lUP + 6il/r + r^u^) + Mr{r - 2M)5'^{r)) ^/3Af Qie(2M - r)5r{r) 



Qci^5r{r) 



V2r9/2 



2r2 






^r^{r-2M) y/r^{r~2M) 
QicSrir) (-6Af2 + SMr + t^lj^) 

2r5 ' '" *"" V ''^ 

Qo {5r{r) (M2 - r'^oj'^) + Mr{r - 2M)5'^{r)) 



^IQ^sl'-^^rir) - f^MQ,f—^SAr) 



(A68) 



(A69) 

(A70) 
(A71) 

(A72) 

(A73) 

(A74) 



2^rii(r-2M) 
where Qx denotes the quadrupole in local frame TSH components. Here Sr{r) is given by the right hand side of (10), 



SJr) = (rci)-^' ,2 ) 
-^ ' V (; 7r3/22'5r(|) 



2 ' ^,y-\ 



(A75) 



18 



d. Quadrupole Components 



Finally, we must obtain the components of the quadrupole from (14). We actually work with the TSH version of 
(14), but this does not pose any problem. For / = 2 the components needed for the present computation read 



Qie = -F{u)2ttV3 Wl 



2M 



M 
2r 



I — Ml 



r-'uj^ 2 + 



3M\ 12M 



Q, = -F(c.)^ (4 + MW) r xMll dr, 



Qe = -F(w)V37r 



" 1- 


. M. ' 

r 


1- 


2M 
r L 



r^w^ ( 2 + 



3M\ 12M 

A 



r J 



X 



-6 



M 



X-[2+'-^\rX' 



5r{r) 



dr, (A76) 



(A77) 



' 2 2 

r u) 



+ 6) -2 



rX' 



Srir) 



dr. (A78) 



The angular integration was already performed. Remember that the Riesz kernel is independent of angular coordinates. 

I 



